Beyond the Jaynes-Cummings model: circuit QED in the ultrastrong couphng regime 
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In cavity quantum electrodynamics (QEDj^^, light-matter interaction is probed at its most fun- 
damental level, where individual atoms are coupled to single photons stored in three-dimensional 
cavities. This unique possibility to experimentally explore the foundations of quantum physics has 
greatly evolved with the advent of circuit QED^ where on-chip superconducting qubits and oscil- 
lators play the roles of two-level atoms and cavities, respectively. In the strong coupling limit, atom 
and cavity can exchange a photon frequently before coherence is lost. This important regime has 
been reached both in cavity and circuit QED, but the design flexibility and engineering potential of 
the latter allowed for increasing the ratio between the atom-cavity coupling rate g and the cavity 
transition frequency cjr above the percent leve P * ^^ ^ ^^ [ While these experiments are well described by 
the renowned Jaynes-Cummings model^^, novel physics is expected when g rea ches a considerable 
fraction of cjr- Promising steps towards this so-called ultrastrong coupling regim j^^^ have recently 
been taken in semiconductor struct ure^^^^. Here, we report on the first experimental realization of 
a superconducting circuit QED system in the ultrastrong coupling limit and present direct evidence 
for the breakdown of the Jaynes-Cummings model. We reach remarkable normalized coupling rates 
g/ur of up to 12% by enhancing the inductive coupling of a flux qubit^^ to a transmission line 
resonator using the nonlinear inductance of a Josephson junctior!^. Our circuit extends the toolbox 
of quantum optics on a chip towards exciting explorations of the ultrastrong interaction between 
light and matter. 



In the strong coupling regime, the atom-cavity cou- 
pling rate g exceeds the dissipation rates and 7 of both, 
cavity and atom, giving rise to coherent light-matter 
oscillations and superposition states. This regime was 
reached in vari ous types of systems operating at different 
energy microwave frequencies, strong 

coupling is feasible due to the enormous engineerability 
of superconducting circuit QED systemJ^. Here, small 
cavity mode volumes and large dipole moments of arti- 
ficial atom^^ enable coupling rates g of aboulP^ 1 % of 
the cavity mode frequency ajr- Nevertheless, as in cav- 
ity QED, the quantum dynamics of these systems follows 
the Jaynes-Cummings model, which describes the coher- 
ent exchange of a single excitation between the atom and 
the cavity mode. Although the Hamiltonian of a realis- 
tic atom-cavity system contains so-called counterrotating 
terms allowing the simultaneous creation ior annihilation 
of an excitation in both atom and cavity mode, these 
terms can be safely neglected for small normalized cou- 
pling rates g/ujr- However, when g becomes a significant 
fraction of cjr, the counterrotating terms are expected to 
manifest, giving rise to exciting effects in QED. 

The ultrastrong coupling regime is difficult to reach in 
traditional quantum optics, but was recently realized in a 
solid-state semiconductor system^^^^. There, quantita- 
tive deviations from the Jaynes-Cummings model have 
been observed, but a direct experimental proof of its 



breakdown by means of an unambiguous feature is still 
missing. In this report, we exploit the potential of fiux- 
biased superconducting quantum circuits to reach the 
ultrastrong coupling regimei^^. For this purpose, we 
increase ^/cjr up to 12% utilizing the large nonlinear in- 
ductance of a Josephson junction (J J) shared between 
a flux qubit and a coplanar waveguide resonator. We 
explicitly make use of the multimode structure of our 
resonator, allowing the direct observation of physics be- 
yond the Jaynes-Cummings model. In equilibrium, our 
transmission spectra reveal anticrossings which can be 
clearly attributed to the counterrotating terms in the sys- 
tem Hamiltonian. These anticrossings are caused by the 
simultaneous creation (annihilation) of two excitations, 
one in the qubit and one in a resonator mode, while an- 
nihilating (creating) only one excitation in a different res- 
onator mode. 

Images of our quantum circuit and a schematic of the 
measurement setup are shown in Figure [ij At a cur- 
rent antinode for the A-mode of a niobium supercon- 
ducting resonator (Fig. [l^-c), a part of the center con- 
ductor is replaced with a narrow aluminum strip inter- 
rupted by a large- area J J. This junction mediates most 
of the inductive coupling between a superconducting flux 
qubit^^ galvanically connected to the strip. The qubit 
consists of three nanometer-scaled JJs interrupting a su- 
perconducting loop, which is threaded by an external 
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FIG. 1: Quantum circuit and experimental setup, a, Optical image of the superconducting A/2 coplanar waveguide 
resonator (light blue box). Black rectangles: area shown in b. Red rectangle: area shown in d. b, SEM image of one of the 
coupling capacitors, c, Sketch of the current distribution of the first three resonator modes. Their resonance frequencies are 
cji/27r = 2.782GHz (A/2, red), uij^i^ = 5.357 GHz (A, blue), and uzj^i^ = 7.777 GHz (3A/2, green), d, SEM image of the 
galvanically coupled flux qubit. The fabrication technology for qubit and resonator is described elsewher^^. The width of 
the center conductor is 20 /xm, that of the constriction 1 /xm. The area of the qubit loop is 180 [ivc? . Orange rectangle: area 
shown in e. Green rectangle: area shown in f. e, SEM image of the large J J with a Josephson inductance Lj, whose large 
inductance is responsible for approximately 85 % of the qubit-resonator coupling, f. One JJ of the qubit loop. The area of this 
junction is 14% of the one shown in e. g. Schematic sketch of the measurement setup. The transmission through the cavity 
at cjrf is measured using a VNA. A second microwave signal at cjs is used for two-tone qubit spectroscopy. The input signal is 
attenuated at various temperature stages and coupled into the resonator (light blue) via the capacitors C^^. The crossed boxes 
represent Josephson junctions. A circulator isolates the sample from the amplifier noise. 



flux bias Scanning electron microscope (SEM) im- 
ages of the qubit loop and the JJs are shown in Fig- 
ure [iji-f- For suitable junction sizes, the qubit poten- 
tial landscape can be reduced to a double-well potential, 
where the two minima correspond to states with clock- 
wise and counter-clockwise persistent currents | ±/p). At 
(^$x = — ^o/2 = 0, these two states are degenerate 
and separated by an energy gap A.^ In the qubit eigenba- 
sis, the qubit Hamiltonian reads ffq = fico'q(Jz/2. Here, 
cjq = Y^A^ + (2/p • 5^^Y /h is the qubit transition fre- 
quency which can be adjusted by an external flux bias. 
We note, that for our flux qubit the two-level approx- 
imation is well justified due to its large anharmonicity. 
The resonator modes are described as harmonic oscilla- 
tors, = huJn{o^t^n + 1/2), where uj^ is the resonance 
frequency and n is the resonator mode index. The op- 
erator dl^ (ttn) creates (annihilates) a photon in the nth 
resonator mode. Due to the inhomogeneous transmis- 
sion line geometry^ (see Fig. [l]i), the higher mode fre- 
quencies of our resonator are not integer multiples of the 
fundamental resonance frequency uji . Throughout this 
work, we refer to the nth mode as the nA/2-mode. Then, 



the Hamiltonian of our our quantum circuit can be writ- 
ten as 

^ = i^q + ^ ^i^n + %n (^1 + «n) (cOS - SlIlO a^) . 

n 

(1) 

Here, ax,z denote Pauli operators, gn is the coupling rate 
of the qubit to the nth cavity mode, and the flux depen- 
dence is encoded in sin 6 = A/hu^ and cosO. The op- 
erator (Jx is conveniently expressed as sum of the qubit 
raising (a+) and lowering (a_) operator. Thus, in con- 
trast to the Jaynes-Cummings model, the Hamiltonian in 
Eq. ([T]) explicitly contains counterrotating terms of the 
form aj^a+ and d^a-. 

Figure [l^ shows a schematic of our measurement setup. 
The quantum circuit is located at the base temperature 
of 15 mK in a dilution refrigerator. We measure the 
amplified resonator transmission using a vector network 
analyzer (VNA). For qubit spectroscopy measurements, 
the system is excited with a second microwave tone cjg 
with power Pg, while using the 3A/2 -mode at uJs/27r = 
7.777 GHz for dispersive readoulPlS^. 

We first present measurements allowing the extraction 
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FIG. 2: Qubit microwave spectroscopy and low power transmission spectra, a, Microwave spectroscopy of the coupled 
qubit-cavity system. The measured transmission magnitude (color coded, blue: low; white: high) is plotted as a function of 
the relative flux bias S^^ and the spectroscopy frequency cJs/Stt. The red broken line indicates the dressed qubit transition 
frequency Cj(p^^^. Inset: center frequency of the qubit spectroscopy signal at = as a function of the probe power Prf • The 
FWHM of the qubit signal is approximately 80 MHz in the low power limit Prf Ps 0. Red line: fit to the linear region^^. The 
green dot indicates the power level at which the spectra in b, c and Fig. [s] are recorded, b, Cavity transmission (3A/2-mode, 
linear scale, arb. units) as a function of and probe frequency c<;rf/27r. c. Same spectrum as in b. Black lines: numerical fit 
of the spectrum of the Hamiltonian ([T]) to the data, d. Cavity transmission (A-mode, linear scale, arb. units) as a function of 
and probe frequency cJrf/STT. The spectrum is recorded at Prf corresponding to n2 ~ 0.9 due to a higher insertion loss of 
this cavity mode. Black lines: numerically evaluated energy level spectrum with parameters from c. 



of the coupling constants of the qubit to the first three 
resonator modes. The spectroscopy data in Fig. [2^ shows 
the dressed qubit transition frequenc}'^ with the ex- 
pected hyperbolic flux-dependence and a minimum at 
(5^x = 0. Furthermore, flux-independent features cor- 
responding to the two lowest resonator modes {ui and 
UJ2) are visible. In principle, a fit to the Hamiltonian 
in Eq. ([T]) would yield all system parameters. However, 
our measurement resolution does not allow us to reliably 
determine the coupling constants in this situation. In- 
stead, we extract from a cavity transmission spectrum 
with negligible photon population. For that purpose, we 
first measure the power-dependent ac-Zeeman shift of the 
qubit transition frequency at S^^ = 0. The data is 
shown in the inset of Fig. [2^. The average photon number 
ns can be estimated using the relation Prf = flshuJsHis^ 
where 1^3/211 ~ 3.7 MHz is the full width at half maxi- 



mum (FWHM) of the cavity resonance and P^f the probe 
power referred to the input of the resonator. Figure (2)3 
shows a color coded transmission spectrum for the 3A/2- 
mode as a function of S^^. The data is recorded at an in- 
put power Prf ^ —140 dBm (green data point in Fig. [2|i, 
inset) corresponding to ns = 0.18. 

We observe a spectrum with a large number of anti- 
crossings resulting from the multimode structure of our 
cavity system. To extract the individual coupling con- 
stants we compute the lowest nine transition frequen- 
cies of the Hamiltonian given in Eq. ([T]) incorporating the 
first three resonator modes. Fitting the results to the 
spectrum of the 3A/2-mode shows excellent agreement 
with the measured data as shown in Fig. We note 
that the spectrum for the A-mode shown in Fig. [2]i can 
be well described without additional fitting using the pa- 
rameters extracted from the 3A/2-mode. For the qubit. 
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FIG. 3: Breakdown of the Jaynes-Cummings model, a, Cavity transmission (3A/2-mode, linear scale, arb. units) 
as a function of and probe frequency cJrf/Svr. Black broken lines in all plots: energy level spectrum obtained from 
the Hamiltonian ([T]). Colored lines in all plots: energy level spectrum obtained from the Jaynes-Cummings model (dark 
blue: l^f, 0,0,1), except for anticrossing region shown in b; green: |e, 1,0,0); light blue: |e, 0,1,0) and |e,2,0,0), which are 
indistinguishable within the resolution of this plot. Grey broken boxes: areas magnified in b-d. b. Single excitation anticrossing. 
The quantitative deviations of the Jaynes-Cummings model from Eq. ([T]) are attributed to a small admixture of the state 
1^,1,1,0). The uncoupled states are indicated by grey broken lines, c. Avoided crossing due to a coupling between the 
degenerate states |^, 0, 0, 1) and |e, 1, 0, 0). This is caused by counterrotating terms in the Hamiltonian (IT]). A detailed analysis 
yields a minor admixture of I ^, 1, 1, 0) (superposition states: ^ l^f, 0, 0, 1)/a/2± (|e, 1, 0, 0)/>/3+ 1, 1, 0)/V6)). This admixture 
has no effect on the reasoning presented in the main text. The energy level spectrum obtained from the Jaynes-Cummings 
model is omitted for clarity, d. Same as c, but for the Jaynes-Cummings model. Within numerical accuracy, no anticrossing 
is predicted, clearly contradicting the data. 



we obtain A/h = 2.25 GHz and 2/p = 630 nA. Most 
importantly, we find coupling rates of gi/27r = 314 MHz, 
g2/27T = 636 MHz, and g3/27T = 568 MHz. The values 
for correspond to normalized coupling rates Qn/^n of 
remarkable 11.2%, 11.8%, and 7.3%, respectively. From 
these numbers, we expect significant deviations of our 
system from a three-mode Jaynes-Cummings model, ul- 
timately proving that we have successfully reached the 
ultrastrong coupling regime. 

In the following, we analyze the features in our data 
which constitute unambiguous evidence for the break- 
down of the rotating-wave approximation inherent to 
the Jaynes-Cummings model. In Figure |3j we compare 
the energy level spectrum of the Hamiltonian in Eq. ([T]) 
to that of a three-mode Jaynes-Cummings model. We 
note that, depending on 5^^^ there are regions where 
our data can be well described by the Jaynes-Cummings 
model, and regions where there are significant devia- 
tions (see Fig. [3^). For our analysis we use the nota- 
tion \q,Ni,N2,N2) = \q) \Ni) (g) iTVs) ^ Ws), where 
q = {g^e} denote the qubit ground or excited state, re- 
spectively, and |A^n) = {lO)? I^)? . . . } represents the 
Fock-state with photon occupation N in the nth res- 
onator mode. At the outermost anticrossings (Fig. [sJd), 
where us ~ cjq, the eigenstates \iIj±) of the coupled sys- 
tem are in good approximation symmetric and antisym- 
metric superpositions of |e, 0,0,0) and |^, 0,0,1). This 
exchange of a single excitation between qubit and res- 



onator is a characteristic of the Jaynes-Cummings model. 
On the contrary, the origin of the anticrossing shown in 
Fig. (Sj) is of different nature: the dominant contribu- 
tions to the eigenstates \iIj±) are approximate symmet- 
ric and antisymmetric superpositions of the degenerate 
states cpi = |e, 1,0,0) and (p2 = |^, 0,0,1). The transi- 
tion from cpi to (p2 can be understood as the annihilation 
of two excitations, one in the A/2-mode and one in the 
qubit, while, simultaneously, creating only one excitation 
in the 3A/2-mode. Such a process can only result from 
counterrotating terms as they are present in the Hamilto- 
nian ([T]), but not within the Jaynes-Cummings approx- 
imation. Here, only eigenstates with an equal number 
of excitations can be coupled. Although counterrotating 
terms in principle exist in any real circuit QED system, 
their effects become prominent only in the ultrastrong 
coupling limit with large normalized couplings gn/^n as 
realized in our system. Hence, the observed anticrossing 
shown in Fig.jS]^ is a direct experimental manifestation of 
physics beyond the rotating-wave approximation in the 
Jaynes-Cummings model. As shown in Fig. [3]i, the lat- 
ter would imply a crossing of the involved energy levels, 
which is not observed. A similar argument applies to 
the innermost anticrossings (see Fig. [3^), although the 
involved eigenstates have a more complicated character. 

In conclusion, we present measurements on a flux- 
based superconducting circuit QED system in the ul- 
trastrong coupling regime. The results are in excellent 
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agreement with theoretical predictions and show clear ev- 
idence for physics beyond the Jaynes-Cummings model. 
Our system can act as an on-chip prototype for unveil- 
ing the physics of ultrastrong light-matter interaction. 
Future explorations may include squeezing, switchable 
ultrastrong couplin^^, causality effects in quantum field 
theor};^^, and the generation of bound states of qubits 
and photons. 
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